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1 Introduction 

Transport, spectral and scattering problems for differential operators on graphs appear frequently in mathe¬ 
matics, natural sciences and engineering [T], [2], |3], |4], [5], |6] and are in the focus of intensive investigations 

mm, 0 , cni, [a, m, mm, m,m- 

In this paper we study the noncompact star-type graph with perturbed radial Schrodinger equation on 
each ray. Due to its physical importance this equation has obtained much attention and we refer for example 
to m, [18] , [19] , [20] , [21] , [22] and the references therein. Our studies dealing with the equation on a metric 
graph can be treated as some continuation of the works [23], [24] (see also the references therein), where 
the important class of the Sturm-Liouville operators having Bessel-type singularity inside the interval was 
investigated. Matching conditions at the vertex arise as a generalization of matching conditions introduced in 
the works above under some additional conditions on behavior of the potential near the singularity. We call 
them the matching conditions of analytical type . One can notice that even for Sturm-Liouville expressions 
with real-valued potentials on finite interval such matching conditions can generate nonself-adjoint operators. 
For this reason we do not restrict our considerations to the case of real-valued potentials and real coefficients 
in matching conditions. 

Our considerations follow in general the scheme recently developed for the investigation of inverse spectral 
problems on quantum graphs. The key point is the solution of ’’partial inverse problem” consisting in 
recovering the potential on one ray from the certain part of scattering data. Although in technical sense 
such problem looks similar to the classical inverse scattering problem for radial Schrodinger equation on 
the half-line, our approach differs from the approach used in the previous works m- Instead of classical 
Marchenko method we use a certain version of contour integral method based in part on the approach 
presented in [25] and in part on the ideas of spectral mappings method [2^. This allows us to provide the 
constructive procedure for solution the inverse problem based a certain linear equation with less restrictive 
requirements on the behavior of the potential at infinity; moreover, contrary to [19] . we are able to prove 
the solvability of this equation. More exactly, under some conditions of regularity and ’’genericity” (which 
are similar to ’’genericity” conditions introduced in [25] and assume the absence of spectral singularities and 
some kind of regularity for small values of spectral parameter) our method works provided the potential 
satisfies at infinity only Marchenko decay condition. 

Main results of the paper are contained in Theorem 5.1 and include the uniqueness result and constructive 
procedure for solution of the inverse scattering problem. 

2 Direct scattering on the semi-axis. Auxiliary facts and nota¬ 
tions. 

Consider the equation: 

~y” + y = = p^y^ ^ ^ (^.i) 
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with complex vq and complex-valued q that will be called as a potential. Let uq = — 1 /4. For definiteness 

we assume that Rez^ > 1/2, i/ and q satisfies the conditions: 

1 OO 

\xq{x)\dx < oo. (2-2) 

0 1 

Here and below we assume := exp(/rlog \z\ + i/iargz), argz € (—7r,7r]. Together with (2.1) we consider 
the unperturbed equation: 

-y" + = \y = p^y, a; > 0 (2.3) 

Denote = 1/2 — z/, /r 2 = 1/2 + Let Cj{x, X), j = 1,2 be the solutions of unperturbed equation (2.3) 
defined as: 

CXD 

Cj{x,X) = x^^ '^CjkX’^x^^, C 10 C 20 = (2z/)"\ 

Cjk = (-l)^CjO I ]^((2s Pj){2.s + Pj - 1) - vq) 

\s=l 

Cj{x, A) are both entire with respect to A. 

Define Sj{x, A), j = 1, 2 as the solutions of (2.1) satisfying the integral equations: 

X 

Sj{x, A) = Cj{x, ~ J (^1 X)dt, 

0 

where g{x, t, A) is a Green function: g{x, t, A) = C'i(a;, A)C' 2 (t, A) — C' 2 (x, A)C'i(t, A). Note that the equations 
are solvable for both j = 1,2 provided the condition (2.2). We also consider the dost solution f{x,p) for 
(2.1) normalized with the asymptotics f^^'>{x,p) = e^P^{{ip)^ + o(l)), x —>■ cx), ^ = 0,1. Such solution exists 
(and unique) under the condition xq(x) G Li(l,oo), analytic with respect to p G D+ := {p : Imp > 0} and 
P~^^ p) is continuous in p G D+. 

The following expansion plays an important role in our further considerations: 

f{x,p) = bi{p)Si{x,X) + & 2 (p)S' 2 (x, A), 

where coefficients hj[p), j = 1,2 are called the Stokes multipliers. The ratio & 2 (p)/&i(p) can be easily 
shown to coincide with the Weyl function m(A), which determines uniquely the potential q. This result and 
constructive procedure for solving the corresponding inverse problem one can find, for instance, in m- 
The following properties of bj{p) arise from the properties of dost solution /(x, p) and the results of [27] . 

Lemma 2.1. For p —>■ 00 the following asymptotics holds: 

bj{p) = P^^ {hf + O(p-i)) , 

where the constants 6°° 0 do not depend on q. 

For p —> 0 the following asymptotics holds: 

b3{p) = P^" {b° + o{l)) 

with some constants bb (depending, in general, on q). 



x^ '^q{x)\dx + 


3 Direct scattering on the graph 

We consider the noncompact star-type metric graph F consisting of finite number of rays emanating 

from their common initial vertex. A function y on the ray TZk we consider as a function of local parameter 
X G [0,oo), where x = 0 corresponds to the vertex. A function y on F we consider as a set of functions 
{yfc}fc=i (where yu = y is considered as a function on [0,oo) as it was specified above). 

On each ray TZk, k = l,p we consider the differential equation 

-y”+{^+dk{x)^y = Xy = p^y. (3.1) 
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We follow the terminology of previous section and use the same notations adding the number k of ray as a 
first index. In particular, we assume Vko = ~ > 1/2, ^ N and qk{x) to satisfy (2.2) (with 

V = Vk)- The complex numbers Vko may be different for different rays but the following restriction is assumed 
throughout the paper. 

Condition v. If Vj ^ Vk then Reiy, ^ Rej^fe. 

Let some function y satisfy (3.1). Then the Wronskians {Ski,y) and (y, Sk 2 ) do not depend on x and we 
can introduce the following linear forms: 


Uki{y) ■■= crk{y,Sk2), 


Uk2{y) ■= ckki{y,Sk2) + crk2{Ski,y)- 

Everywhere in the sequel we assume ak ^ 0, <7k2 7 ^ 0, fc = l,p. 

Now let y = {ykYk=i ^ function on T such that each of yk satisfies (3.1). We define the matching 
conditions at the vertex of T as follows: 

p 

Ujiiyj) = Uki{yk),j 7 ^ k, ^Uj 2 {yj) = 0. (3.2) 

1=1 


Matching conditions (3.2) are generalization of matching conditions introduced in |23] for Sturm-Liouville 
operators having singularities inside the interval. They also can be considered as some generalization of 
standard (i.e. continuity together with Kirchhoff condition) matching conditions for quantum graphs. 

With each ray TZk we associate the Weyl-type solution ipk (p) that we determine as a function on T tpk (p) = 
{ipkjix, such that: 

• each of ipkj, j = 1,P solves the differential equation (3.1); 

• ipkj{x, p) = O as a: —7 oo for j ^ k; 

• tpkkix, p) = e“*^“(l + o(l)) as a: oo; 

• i/’fc(p) satisfies the matching conditions (3.2). 

Let k be arbitrary fixed. We look for 'tpk in the following form: 


i’kj{x,p) 


lkj{p)fj{x,p),j ^ k, 'tpkkix,p) = Jkk{p)fk{x,p) + 


2ip 

bki{p) 


Sk 2 {x,X). 


(3.3) 


It is clear that any function of the form (3.3) with arbitrary coefficients "fkj satishes all the conditions 
determining the Weyl-type solution except matching conditions (3.2). Substituting (3.3) into (3.2) and 
taking into account that (/j, Sj 2 ) = bji, {Sji, fj) = bj 2 we arrive at the following linear algebraic system: 


p 

(Tjbji-ikj + I3k=0, j = l,p, '^{cTjibji + aj2bj2)lko 

i=i 


—<Xk25k, dk 


2ip 

bki 


(3.4) 


with respect to the values ^kj^j = 1,P, Pk- Cramer rule being applied to (3.4) together with representation 
(3.3) shows that ij^kip) exists (and is unique) for all p G SI 4 . outside the zeros of the function bki{p)A{p), 
where 

p 

^{p) = ^(o'si&si(p) + crs2bs2{p)) Y[<Jjbji{p) (3.5) 

S = 1 j^s 

is a determinant of the system (3.4). In the sequel we assume that the following ’’genericity” condition holds. 

Condition G. All the functions bki(p)A(p), k = l,p have no real zeros (with possible exception of 
p = 0 ) and no multiple zeros. 

m 

Let {I, ...,p} = (J where the pairwise disjoint subsets are such that for any j, k G one has 
«=i 

Pj = k'k = and Reri > .. .Rer^. The following properties of characteristic function A(p) are the direct 
sequences of representation (3.5) and lemma 2.1. 
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Lemma 3.1. For p —>• cx) the following asymptotics holds: 


C=i 


where 

p 

af = Y^ fJs2b‘^Y[ajbJ^, pi 

sG/j j^s j=l 

For p —> 0 one has: 

A(p) = p'"i(d° + o(l)), (f = ^(cr^i 6 °i + cj, 2 b° 2 ) n 

S = 1 j^s 

It is important to notice that the coefficients a|° depend on the numbers Vk, k = l,p and coefficients 
CTfc, CTfci, CTfc 2 in matching conditions (3.2) and do not depend on the potentials qu- In what follows we assume 
that the following conditions of ’’regularity” and ’’genericity at 0 ” hold. 

Condition i?oo- ^ 0. 

Condition Rq. bj^ 0, j = 1,P, 7 ^ 0. 

Our further considerations deal with ipkk{x,p). Using (3.3) and the representations: 

Ikkip) = Skip) • n (3-6) 

that can be obtained from (3.4) using the Cramer rule we deduce the following assertion. 

Lemma 3.2. Under the Condition Rq one has 

Ikkip) = O 

as p —>■ 0. Moreover, 

i’kk^\x^p) = o (p^""), C = 0,1 

as p —5- 0 for any fixed a; > 0. 

It follows from lemmas 2.1, 3.1 that under the conditions G, Rq and Roo the function 6fci(p)A(p) has at 
most finite set of simple zeros, ipkkixjp) has at these points either a simple pole or a removable singularity. 
Denote the set of poles of tpkkixjp) in as . It follows from representation (3.3) that for any po € Z'^ 
one has: 

resp^pgipkkix, p) = akipo)fkix,po) = afc(po) exp(ipox)(l +o(l)),x cx). (3.7) 

with some constant akipo)- Indeed, (3.7) is obvious if po is a zero of A(p); if po is a zero of &fci(p) we take 
into account that Sk 2 ix,po) is proportional to /^(x, po) and arrive again to (3.7) (let us recall that under 
the condition G A(p) and bj-iip) have no common zeros). 

Further, for real nonzero p usual arguments yield the following representation: 

'ipkkix,p) = fkix,-p) + rfc(p)/fc(x,p) = exp(-ipx) + rfc(p) exp(ipx) + o(l),x -7 00 . (3.8) 

We call rfc the reflection eoejflcient associated with the ray TZk, the set Jk '■= {rt, Z'^, akip), p G Z ^} is called 
the scattering data assoeiated with the ray TZk, the set J = {Jk}k=\ is called the scattering data for F. 

4 Inverse scattering. Partial problem. 

This section is devoted to the following partial inverse problem. 

Problem /P(fc). Given Jk, fc G {I,... ,p} find Qkix), x > 0. 
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Let L denotes the problem on L consisting of differential equations (3.1) on each ray TZj, j = l,p and 
matching conditions (3.2). Together with L we consider the problem L with equations of the same form 
(3.1) and matching conditions (3.2) but different potentials qj, j = l,p. We assume that coefficients in (3.2) 
are the same for L and L and Moreover, we assume that the conditions G, Rq and i?oo hold 

for both L and L. We agree that if ry denotes some object related to L then fj denotes an analogous object 
related to L and rj := ij — fj. 

In what follows we shall also use the following notations. If A denotes some matrix then Aj denotes 
its j-th row. If / is some function holomorphic in deleted neighborhood of po then f{m)iPo) denote the 
coefficients in the Laurent expansion: 

OO 

f{p)= {P- PoTf{m){Po)- 

m— — oo 

If f(p) is some function meromorphic outside the real axis then for real p we denote f±(p) '■= lim f(p±is) 

£->+0 

From this point in this section we assume that k is arbitrary fixed. The key role in our further consider¬ 
ations is played by the spectral mappings matrix P{x,p) := '^{x, p), where: 

'i!{x,p):= /^(xp) > d e ^'(x,p) := ^'(x,-p), p e := {p : Imp < 0}. (4.1) 

Lemma 4.1. For each fixed x > 0 P{x, p) is bounded for p -5- 0 and p —>■ oo. Moreover, for p —>■ oo one 
has Pi{x,p) = Ii + 0{p~^), where (and everywhere below) I denotes the identity matrix. 

Proof. By virtue of symmetry it is sufficient to consider p from upper half-plane 12+. For nonzero 
p G 12+ one has det = det ^ = 2ip and: 

2ipP^i{x,p) = 'iljlY\x,p)fl,{x,p) - fk~^\x,p)iJkk{x,p), 

2ipP^2{x,p) = fjf~^\x, p)ipkkix, p) - ■ip[Y\x, p)fkix, p). 

Lemma 3.2 together with the estimate fY^\x, p) = O (|p^'“i |) provide the boundedness of P{x, p) for p —> 0. 
For p —>■ oo we use the asymptotics: 

fk^\x,p) = (ip)^ exp(ipa;)[l], (a;,A) = l3koP~^'^'" ((-ip)^ exp(-ipa;)[l] -f- (ip)^ 7 fco exp(ipa;)[l]) , (4.2) 

where [1] := 1-1- 0{p~^) and the constants /3ko, 7fco depend only on i^ko (and therefore are the same for L 
and L). In particular this yields the estimates: 

(x, p) = O {iip)^~^ exp{ipx)) , (x, X) = O exp(-ipx)) . (4.3) 

Moreover, lemma 3.1 provides the estimates: 

A(p) = O (p'^i+2^1-1) , A-i(p) = O . (4.4) 

On the other hand one has the asymptotics (following directly from (3.5)): 

Afc(p) = p^-'^-(dr + 0(P-')), dT=alX{ a,bf,. 

Together with lemma 2.1 and the estimates (4.4) this yields: 

7fcfc(p) = O (p2i2M-2Ti^ ^ ^ ^ 0{p-^). (4.5) 

The asymptotics (4.2), (4.3), (4.5) provide the estimates: 

V'lfc (x, p) = 0 (p^ exp(-2pa;)) , (x, p) = O (p«"^ exp(-ipa;)) . (4.6) 

Together with (4.2), (4.3) the estimates (4.6) yield the required boundedness of P{x,p) for p —>■ oo and the 
estimate Pi 2 {x,p) = 0{p~^). In order to obtain the estimate for Pii(a;,p) — 1 we write: 

2ipPii{x,p) = il)kk{x,p)fk{x,p) - fk{x,p)tpkk{x,p) = 2ip + ipkki.x, p)fk{x, p) - fk{x,p)'tpkk{x,p) 
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(where we take into account that {tpkk,fk) = ‘^ip) and use again (4.2), (4.3), (4.6). □ 

It is clear that P{x, p) is meromorphic in p outside the real axis and (under condition G) for all nonzero 
real p has the limits P±{x, p) which are continuous and bounded on ]R\ {0}. All the possible poles of P{x, p) 
are necessarily belong to the set Zk := Zfc U Zfc, where Zk = {±/o, p € Z^ }. 

Lemma 4.2. Any po € Zk is either a simple pole or a removable singularity for P{x,p). Moreover, the 
following representation holds: 

P(_i>(a;,po) = ^'{o>(a;,Po)^(do)(^'"^)(o)(a:,Po), 


where v{po) := 0 ii po ^ Zk and 


v{po) 


0 ak{po) 
0 0 


if po € Zk- Here ak{po), po & Z'^ are the constants from (3.7) and ak{po) '■= —ak{—po) for po G Zk H H_. 

Proof. It follows from (3.7) that po G Zk is (at most) a simple pole for 4'(a;, p) and the following relation 
holds: 

4'(_i)(x,po) = ^'(o)(a:,Po)v(do)- (4.7) 

Moreover, since det = det ^ = ±2ip for ±p G H+, the matrix ^~^(x, p) has (at most) a simple pole in po 
and from vj/vp-i = J one can easily obtain: 


^)(-i>(x,po) =-v(do)(^ ^)(o>(a;,po)- 


(4.8) 


Thus, P{x,p) has in po a pole of multiplicity 2 or less. But P{- 2 ) = ^(-i)= ~'f'(o>'^'5(^~^)(o>• 
We notice that because of special structure of matrices v{po) and v{po) we have v{po)v{po) = 0 and therefore 
P{- 2 )ix, Po) =0. Now we calculate P^_i) = + ^(o)(^~^)(-i> and using (4.7) and (4.8) obtain 

the required representation. □ 

Our first result in this section is the following uniqueness theorem. 

Theorem 4.1. Jk = Jk implies qk{x) = qk{x) for a.e. a; > 0, i.e. specification of the scattering data 
associated with some ray TZk uniquely determines the potential qk on this ray. 

Proof. First we notice that (3.8) yields the following relation for the matrices 'I'±(a;,p): 


'i'+{x,p) = 'i'-{x,p)v{p), peIR\{ 0 } 


(4.9) 


with ’’jump matrix” 


v{p) 


rk{p) I 

I - rk{p)rk{-p) -rk{-p) 


(4.10) 


From (4.10) one can notice that fk = Xk implies v{p) = v{p) for all p € R \ {0} and therefore P+(x,p) = 
P_(x,p) for each fixed x > 0 and all nonzero real p. Thus, if Jk = Jk the matrix P(x,p) is meromorphic 
in C \ {0} with possible poles in points of the set Zk- Moreover, Jk = Jk implies Zk = Zk = Zk and 
f’(po) = v{po) for any po G Zk- Then, by virtue of lemma 4.2 we have P(_i)(x,po) = 0, i.e. P(x,p) has 
actually removable singularities in all po G Zk- Thus, P(x,p) is holomorphic in C \ {0} and by virtue of 
lemma 4.1 bounded for p —> 0 and p —> oo. Moreover, by virtue of lemma 4.1 Pi(x,p) — Ii vanishes at 
infinity. This means that Pi(x,p) — Ii is actually equal to 0 identically. Thus, we have Pi(x,p) = Ii that 
yields fk{x,p) = fk{x,p) and subsequently qkix) = qk(x) for a.e. x > 0. □ 


Our next goal is a constructive procedure for solving the problem IP{k). From this point we assume that 
the numbers Vjo^ j = l,_p and coefficients of the forms Uji, Uj 2 in matching conditions (3.2) are known. Also 
we suppose that L is chosen a priori ’’model” (or ’’reference”) problem with known potentials qj,j = l,p. 
Below we reduce the problem IP{k) to a certain linear equation (’’main equation”) in some Banach space 
and prove the unique solvability of this equation. 

We consider again the spectral mappings matrix P(x,p). By virtue of lemma 4.1 one has /|^|^^(/^ “ 
p)~^dp{Pi{x, p) — Ii) —5- 0 as P —5- oo that yields the following basic relation: 


OO 

Pl(x,p)-/l= ^ (p-p)”^Pl,(_l)(x,p) + (Pl,+ (x,p) - Pi _(x,p)) (4.11) 

, 27TI J — p 




6 



with arbitrary p S C \ (M U 

Let us introduce the matrix V{x,p) := PZ^{x, p)P+{x, p) = ^-{x,p)v{p)^^^{x,p), p G R \ {0}. We 
notice that V is nniquely determined by L and scattering data Jk- We also note that in view of lemma 4.1 
and det P = 1, V{x,p) is continuous and bounded with respect to p € R \ {0} and any fixed a; > 0. 

We define 

di{x,p,po)= [(p-p)”^4'“^(a;,p)] , . , p GC\Zk,po G Zk 

d2{x,po,p)= [(p-p)"^4'(a;,p)] , . , p G C\ Zk, po G Zk 

D(x,p,p) := {p- p)~^{'i>~^){x,p)'i>{x,p), p^ p,p,pG<C\Zk, 

D{x,p,po) := [D{x,p,p)],qs , p GC\Zk,po G Zk 

D{x,Pq,p)-.= [D{x,Pq,p)],q. , p GC\Zk,po G Zk 

' ' P=P0 

Dix,po,po) ■■= [D{x,p,po)],q) = [D{x,po,p)],ox , Po € Zk, po G Zk 

P=P0 P=P0 

and dj{x,p,p), j = 1,2, D{x,p,p) by similar formulae, where 'I'(x,p) is replaced with 4'(x,p). We also 
define: 

A{x,p,p) := D(x,p,p)v(p), A(x,p,p) := -D{x, p, p)v{p), p G Zk, pG Zk- 
Lemma 4.3. For p G Zk, p ^ Zk one has: 

{P - P{-i){x, p) = 4'(o>(a;,p)t)(p)(Ii(x,p,p), 

{p-p)~^P(_i){x,p)^{x,p) = 'H(Q){x, p)v{p)b{x, p, p). 

For ^ G Zk, p,p ^ Zk 

ip - - P)~^'^~^ix, p)P{-i) [x, p) = D{x, p)v{0Dix, p, 0- 

Proof. The relations can be easily verified by direct calculation. Consider, for instance, 

di{x,p,po)= {p- p)~^b~^{x,p) = 

^ J P=Po 

(p-po)"^(^"^){o>(a;,Aio) + [ip-p)~^](i) ■ i^~^){-i)ix,Po)- 

' ' P=P0 

Multiplying this by ^q){x, po)vipo), then using lemma 4.2, relation (4.8) and taking into account that 
v{po)v{po) = 0 we obtain the first of the required relations. The others can be obtained in a similar way.□ 

The assertion below follows directly from definition of dj{x, p, p). 

Lemma 4.4. For any fixed p G Zk di(x, p, p) is holomorphic in p G C \ In particular for real p one 
has di{x, p -\- iO, p) = di{x, p — iO, p) = di{x, p, p), moreover, di{x,-,p) G L2 (R,C^). Analogously, for any 
fixed p G Zk d2{x,p,p) is holomorphic in p € C \ Zk, d2{x,p, •) G L2 (R, C^). Here and below elements of 
are considered as row-vectors. 

In our further considerations we will also use the Cauchy operators: 

00 

iCf){p)-.= ,^ f -^fip), pGC\R, 

2 tti J p- p 

— 00 

iC±f)ip) := (C/)±(p), pgr 

and the Plemelj-Sokhotskii formula C± = P ± E/2. Here (and everywhere below) E denotes the identical 
operator and 

00 

iPf)ip) = ,^ f -^fip), peR 

2 tti J p- p 
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where a principal value of integral is assumed. We recall that C± and V are bounded operators in ^ 2 ( 8 ) M) 
(for any finite-dimensional space M) and the following relation holds for two arbitrary matrices Fi,F 2 of 
suitable orders with quadratically integrable elements: 


00 

J dfi{C±Fi){^i)F2ifi) 

— 00 


00 

- J dnFi{fi){C^F2){n). 

— 00 


Define $(x,p) := Pi^+{x,p) — Pi-{x,p) for real p and $(x,p) := (a:,p)fi(p) for p G Zk- One can 

easily notice that by virtue of lemma 4.1 •)|g G L 2 (K., C^). Now we return to the relation (4.11) and 

using the notations introduced above and lemma 4.3 rewrite it as follows: 

Pi{x,p) - h= ^ <^{x,p)di{x,p,p) + {C‘^{x,-)){p), / 9 GC\(IRUZfc). (4.12) 

Taking in (4.12) the limits as ilmp - 5 - 0 and substituting them into the jump relation — Pi.-V = 0 we 
arrive at: 

00 

^(a;,p) Q 4 +^l^(a:,p)^ -f ^ J p){I - V{x, p))+ 

— 00 

^ <P{x,p)di{x,p,p){I-V{x,p)) + h-Vi{x,p) =0, pGlR\{0}, (4.13) 

fiGZk 

where a principal value of integral is assumed. In the particular case when = 0 the relation (4.13) can be 
considered for each fixed a: > 0 as a linear equation with respect to $(x, •). In general case (4.13) should be 
completed with some relations at the points p G Zk- We proceed as follows. Multiplying (4.12) by ^{x,p) 
and using lemma 4.3 we obtain: 

'^i{x,p) -^i{x,p) = ^ p)v{p)D{x, p, p) + {C<^{x)) {p)^{x, p), pGC\(]RUZfc). 

fiGZk 


Now we take an arbitrary po G Zk and multiply the relation above by v(po)- Thus we get: 

•i>i(x,p)v(po) - ^i(x,p)v(po) = 

00 

'^i,(o}(x,p)v(p)D(x,p,p)v(po) + ^ J p)^{x, p)v{po). (4.14) 

l-t'^Zk —00 

Consider again the relation (4.14). Taking the coefficient [. ■ ^(o) \p-p Laurent series of its both sides and 
using lemma 4.3 we arrive at: 


*^{x,po) 


2Tri 


dp^{x,p)d 2 {x,po,p)v{po) - ^ ^{x,p)A{x,po,p) = '^i^i^o)ix,po)v{po),po e Zk, (4.15) 

fi^iZk 


where po G Zk is arbitrary. Together with (4.13) (4.15) forms the complete linear system with respect to 
^{x,p), p € M.L) Zk- We can summarize our considerations as follows. 

Theorem 4.2. Let dj{x,p,p) and A{x,p,p) constructed as described above by given scattering data Jk 
and given model problem L. Then 


<^{x,p) 


Pi,+ {x,p)-Pi-{x,p), pGM, 
'^i,{o){x,p)v{p),p € Zk 


for each fixed a; > 0 solves the linear system (4.13), (4.15). 

In order to complete our procedure of solving the problem IP{k) we are to show the unique solvability of 
the specified linear system. First we rewrite the system as a linear equation in Banach space "H := T-lr (BPa, 



where 'Hr = -^2 (R, C^), 'Hd = As we already mentioned above for each fixed x > 0 •) can be 

considered as an element of 'H. Thus, the system (4.13), (4.15) can be written in the form A4> = G, where 


G{p) := 


yi{x,p) - ii, p e K, 
^i,{0)ix^p)v{p),p G Zk, 


while ,4. is a linear bounded operator in 'H defined by the following operator matrix: 


A = 


A rr 

>^dr >^dd 


Arr'P = C+ip-{C_ip)V, {ArdP){p)= 'P{p)dl{p, p){I - V {p)), p 


fd'GZk 


{Adr<p){p) = j p{0d2{p,0vip)d^, p€Zk, iAddp)ip) = p{p) - ^ p{p)A{x,p,p), p G Zk- 

— 00 

Here and below we assume that parameter x > 0 is (arbitrary) fixed and for the sake of brevity we omit 
it in all the argument’s lists. Boundedness of the operator A follows from lemma 4.4 and boundedness of 
operators C±, function V. Now we can provide the following improvement of Theorem 4.2. 

Theorem 4.3. For each fixed x > 0 4’(x, •) is a unique solution in 'H of the equation A(t) = G. Moreover, 
operator A has a bounded inverse operator A~^ = B, where: 


B = 


Brr Brd 

Bdr Bdd 


Brrf = C+{fP+)P+-G-{fP+)P-, {Brdf){p)= ^ / (p)dl (p, p) (P+(p) - P-(p)), p 


l-i'GZf- 


00 

(Hdr/)(p) = -^ J f{ 0 P+{ 0 d 2 ip, 0 v{p)d^, pGZk, {Bddf){p) = f{p) - ^ /(/i)A(x,p,p), p G Zk- 

— 00 f-t^Zk 

Here and below P := P~^ while P denotes as above the spectral mappings matrix. 

Lemma 4.5. The following relation holds for nonreal p ^ p, p, p G C \ Zk'- 


1 


d^ 


-{p+io-p-io) = -^p{p) - -^pip) + E d2{^,p)mMp,o- 

p-p p-p 


27ri J (p-Oi^-p) 

— 00 

For p € C \ (R U Zk), p G Zk one has: 

00 

^ f -^di{^,p)iP+y - P-iO) = di{p,p) - di{p,p)Pip) - ^ A{^,p)di{p,^). 

2tti J p-^ y- 

For p G C \ (R U Zk),p G Zk one has: 

00 

f T^iP+y - P-y)d 2 {p,y{p) = P{p)d 2 {p,p)v{p) - d 2 {p,p)v{p) + E d2i^,p)v{0Mp,0- 

2tti J ^-p y- 

For p,p G Zk one has: 

00 

^ J dy,p){Py) - Py))d2{p,C)v{p)dS,= A{p,p) + A{p,p) - ^ A{^,p)A{p,^). 


ieZk 
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Symmetrically, one can obtain the following relations: 




— -iP+iO-P-iO) = - -^Pip) - E d2{pp)H0diip,0 

— //I n — II n — II • ^ 




2Tri J (p-Oi^-p) P-P P-P 

— OO 

for p, /X G C \ (R U Zk), 

OO 

7^ f -^di{^,p)iP+iO - P-iO) = diip,p) - diip,p)Pip) - ^ A{^,p)di{p,^) 

2m J p-^ 

— OO qGZk 

for p € C \ (R U Zk), p G Zk, 

OO 

f j^{P+iO - P-{0)d2ip,0v{p) = P{p)d 2 {p, p)v{p) - d 2 {p,p)v{p) + E d2iPp)v{0Aip,0 

2m J ^-p ^ 

for p G C \ (R U Zk),p G Zk, 

OO 

/ dAPp){P^{P)-P-{P))d 2 {p,i)v{p)di = A{p,p)+A{p,p)- A{^,p)A{p,^) 




for p, p G Zk. 

Proof. All the relations can be obtained in a similar way based on the relation: 

di 


J (p-e)(C-p) 

\i\=R 


P(e) = 0, p,pGC\(RUZfe), 


(4.16) 


that follows from lemma 4.1. 

Let us show how to get, for instance, the forth relation. Multiplying (4.16) by 4'“^(p) from the left, by 
4' (p) from the right and applying the residue theorem we obtain: 


d^ 


1 


27ri J (p-Oi^-p) 


^-\p)iP+iO - P-{0)^ip) = 


E 7--r«'-'(M)^^(-i>(e)^(p) + ^^'"'(A^)^'(p)-^'I'"'(A^)'I'(p), p,pGC\(RUZfe). 

(p-Oi^-p) p-p p-p 


Using lemma 4.3 we rewrite this as follows: 


d^ 


1 


27ri J (p-Oi^-p) 


vI/^i(p)(P+(e)-P_(0)4'(p) = 


E p)v{OD{p, 0 + D{p, p) - D{p, p), p, p G C \ (R U Zk). 

i&Zk 

Taking (for arbitrary fixed p G C \ (R U Z^,), po G Zk) the coefficients [...](o) I _ in Laurent series of both 
sides we arrive at: 

OO 

i / ■^^d^i^,Po)ip+io-p-imiP) = 

— OO 

- E Po)PiP, 0 + P{p^ Po) - D{p, Po), Po G Zfe, p G C \ (R U Zk). 

i&Zk 
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Finally, taking (for arbitrary po,l^o € Zi~) the coefficients _ in Laurent series of both sides and 

multiplying them by v{po) we obtain the required relation. □ 


Proof of Theorem 4.3. Boundedness of the operator B follows from lemma 4.4 and boundedness of 
operators C ±, functions V, P±,P±. The proof consists in the direct calculations of elements of the operator 
matrices AB and BA. 

Position (r,r). Let ip = Brrf- We rewrite it as ip{p) = f{p) + C-{fP+) (p)(P+(p) — P_(p)). Then 
{Cip){p) for arbitrary fixed nonreal p G C \ Zk can be written in the following form: 


iC<p){p) 


iCf){p) - 


2Tri 


C-{fP+) {p)Fp{p)dp,, 


where Fp{^) := {p — p) ^(P+(/x) — P_(//)). Using the relation 


oo oo 

-/ iC-Fi){p)F2{p)dp.= J Fi{p)(C+F2){p)dp 


we rewrite it as follows: 


{Cip){p) = (C/)(p) + — / ifP+)ip)iC+Fp)ip)dp. 


Now let us consider {CFp){p). For nonreal p GC\Zk using lemma 4.5 we calculate: 


iCFp)ip) = 


d^ 


2-Ki J {p-Oii-d) 


— -AP+io-p-iO) = -^pid) - -^pip)+ E d2{^,pm)MpA)- 

— //I n — II n — II f ^ 


p-p 


p-p 


CeZk 


Taking the limit as Im/i —>■ +0 (while p S C \ (K U remains fixed) we obtain: 

{c+Fp)ip) = ^p+(p) - -^Pip) + E d2i^,p)mMp,o- 

n — II n — II • ^ 


p-p 


p-p 




Thus, we can write: 


iCip){p) = (C'/)(p) + :F f (/P+)(p) i -^—p^{p) - —P{p) + E d2i^iP)v{0diip,0 \ dp = 

2Tn J \ P- P P- P f-A 

— OO 1. ) 


— j ifP+)ip) I E d2{^,pm)MpA) > dp+ [cifP+)\ {p)P{p) 

— oo J 

that yields the following relation: 

(ArrBrrfKp) = (ArrPKp) = \C+{fP+)] {p)P+{p) - \C.ifP+)] (p)P_ (p)T(p) + 


F J (/p+)(p) I E d2(e,M)*(e)rfi(p,o| {I-Vip))dp. 




On the other hand we have: 


(ArdBdrfKp) ^ E {-2^^ / f(.P)P+ip)d 2 i^,P)mdp \ diipAKl -V{p)) 
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and thus arrive at: 

iArrBrrf){p) + {ArdBdrfKp) = [C+{fP+)] {p)P+{p) - [C-{fP+)\ {p)P.{p)V {p). 

Taking into account that P+ = P-V and P = P~^ we obtain finally {ArrBrrf){p) + {ArdBdrf){p) = 
[{C+-C.){fP+)]{p)P+{p) = f{p). 

Now let / = ArrP- Then {fP+){p) = {(pP+){p) + {C-ip){p){P+{p) — P-{p)). Consider the function 
[C{fP+)]{p). Proceeding as above we obtain: 

OO 

[C{fp+)\{p) = [C{pp+)\{p) + ^ / p{p)[C+Fp\{p)dp, 

— OO 

where Pp{p) := {p — p)~^{P+{p) — P-{p)). Using lemma 4.5 we calculate: 

[C+Fp\{p) = -^P+{p) - -^P{p) - Y. d2{^,Pm)diipA), 


[C{fP+)]{p)P{p) = {Cip){p) - ^ y < Y d2ipp)v{0diip,0 } P{p)dp 

— OO J 

that yields finally: 

{BrrArrP)(.p) = {Brrf){p) = p{p) “ ^ < Yj d2iP P)viOdl{p, 0 > (P+ip) “ P-{p))dp- 

^*_oo J 

On the other hand we have: 

{BrdAdrP){p) = Yj |^/ ‘fi(B)d2{^, p)v{Odp\ di{p, ^){P+{p) - P-{p)) 


and thus {BrrArrP){p) + {BrdAdrP){p) = Pip)- 

Position (r,d). Calculate ArrBrd + ArdBdd- Let p = Brdf- Then 

Arrp{p) = Y [C+Rp - {C-Rp)V] {p), 


where Rp.{p) '■= di{p, p){P+{p) — P-{p)). Consider for nonreal p £C\Zk the function: 

OO 

— OO 

Using lemma 4.5 we rewrite this as follows: 


CRpip) = -di{p,p)+di{p,p)P{p)+ Y M^^P)di{P,0- 


Taking the limits ±Imp —>■ 0 we obtain: 


[C+Rp - {C-Rfj,)V] ip) = -diip,p)il -Vip))+diip,p)iP+ip) - P-{p)V{p)) 

+ Y Ai^,p)d,ipA)iI -Vip)). 


Thus, taking into account that P+ = P_U we arrive at: 


iArrBrd)ip) = { Y f(B)FiP,P) \ (I ” Vip)) 
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where 


Fip,fJ.)= ^ A{^,p)di{p,^) - di{p,p,). 

i&Zk 

On the other hand direct calculation yields: 


{ArdBddf){p) 



fiO- 




{i-v{p)), 


that can be rewritten as 

iArdBdd){p) = - I E fiA)F{p,p) \ {I-V{p)). 

Thus, we have ArrBrd + ArdBdd = 0. Symmetrical calculations show that BrrArd + BrdAdd = 0. 

Position (d,r). Calculate AdrBrr + AddBdr- Let p = Bdrf- Then Addp{p) can be written in the 
following form: 


Addp{p) 


OO 

i / fMPM 

— OO 


-d2{p,p)vip) + E d2i^^p)v{OMpA) 


dp. 


Consider the function Rp{p) := {P+{p) — P-{p))d 2 {p, p)v{p) (where p G Zk considered as a parameter). By 
virtue of lemma 4.5 one has for nonreal /i G C \ (R U Z^): 

(CRp){p) = P{p)d 2 {p,p)v{p) - d 2 {p,p)v{p) + E d2{£,,p)v{0Mp,0- 

i&Zk 

Taking the limit as Imp —>■ +0 we obtain for real p: 

{C+Rp){p) = P+{p)d 2 ip,p)v{p) - d 2 {p,p)vip) + ^ d2i^,p)v{0MpA)- 

i&Zk 

Thus, we can write: 


AddPip) = ^ / fit^)P+id) [iC+Rp)ip) - P+{p)d2ip,p)v{p)^ dp = 


TT- / fid)P+{p)iC+Rp)ip)dp- — / fip)d 2 ip,p)v{p)dp. 


2Tri 


2TTi 


Using the relation 


I OO 

Fi{p)(C+F2){p)dp = - J {C_Fi){p)F2{p)dp 


we rewrite this as follows: 


{AddBdr f)p = Addp{p) = — 


27ri 


C_(/P+) {p)Rp{p)dp- 


2711 


f{p)d 2 {p,p)v{p)dp = 


1 

27ri 


C-{fP+) {p){P+{p) - P-{p))d 2 {p,p)v{p)dp-— / f{p)d 2 {p,p)v{p)dp. 


Now let p = Brrf- We rewrite it as p{p) = f{p) + C_ j^(/P_(.)J {p){P+{p) — P_(/x)). Then 

{AdrBrrf^ip} — {AdrP^{p') ~ 
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^ I fi^^)d2ip,^J')v{p)d^i+I \C-{fP+)\ {fi){P+{fi) - P-{fi))d2ip,p)v{p)dn. 


Analogous calculations yield 


{BddAdrP)ip) = ^. f [C-ip]ip)iP+{p) - P-{p))d 2 ip,p)v{p)dp +^ I ip{p)P+{p)d 2 {p,p)v{p)dp. 


On the other hand, if / = ArrP then {fP^){p) = (p{p)P+{p,) + [C-(p]{p,){P+{p,) — P-{p)) and 

{BdrArrP)ip) = {Bdrf){p) = 

oo oo 

/ 'd^^B)P+{p)d 2 {p,p)v{p)dp,- j [C-(p]{p){P+{p) - P-{p))d 2 {p,p)v{p)dp. 

— OO —OO 

Thus, we have BddAdr + BdrArr = 0. 

Position (d,d). Let p = Brdf ■ Then 

{Adrp){p) = ^. J I ^ f{p)di{^,p) \ {P+{^) - P-{^))d2ipA)v{p)d^- 

Using lemma 4.5 we can rewrite it as follows 


{AdrBrdf)ip) = iAdrP){p) = ^ /( m ) < Mp, P) + A{p, p) - ^ A{^,p)A{pA) \ ■ 


On the other hand we have: 


{AddBddf){p) = fip) - ^ f(p)A(p,p) - ^ /($) - f(p)A{i,p) \ A(pA) 

^GZk \ fd-^Zk J 

and thus we obtain {AdrBrdf){p) + {AddBddf){p) = f{p)- Symmetrical calculations yield BdrArd + BddAdd = 

E. □ 

Now we can formulate the constructive procedure for solution of the problem IP{k). 

Procedure 4.1. Given scattering data Jk (with some fixed k & {1,... ,p}) and voj, Uj,aji , (7j2, j = I,p. 

1. Choose L with the same VQj,aj,aji,(jj 2 , j = l,p satisfying the conditions G, i?o, Roo- 

2. Calculate v{p), p G K U using (4.10) and lemma 4.2. 

3. From given u(p), p G R U and L find V(x, p), dj(x, p, p),j = 1,2, A(x, p, p). 

4. For each fixed a; > 0 find 4>(x, p), p G M U as a (unique) solution of the linear system (4.13), (4.15). 

5. Given $(a;, p) calculate Pi(a;, p), a; > 0, p G C \ R U via (4.12). 

6 . Given Pi (a;, p) find/fc(x,p) = Pii{x, p)fk{x, p)+Pii{x, p)jl{x, p), ipkkix,p) = Pii(x, p)^fcfc+^'ii(a;, p)'i/'fcfc(a;, p). 

7. Calculate qkix) = f'^{x, p)/fk(x, p) + p^ — (where p for each fixed x > 0 is arbitrary such that 

/fc(x,p) ^ 0). 

5 Inverse scattering on the graph 

Here we consider the following ’’complete” inverse scattering problem. 

Problem /P(F). Given scattering data J (= {^fcl^Ii) hud the potential on F, i.e. all the functions 
qk, k = l,p. 

First we note that the procedure described in previous section allows us to recover (uniquely) the poten¬ 
tials < 7 fc, /c = l,p — 1. Thus, in order to complete our solution of the problem we have to recover qp. Let us 
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consider the matching conditions for Weyl-type solution ipi{p). The following relation is a direct sequence 
of (3.2): 




(5.1) 


Since for j = 2,p ipkjix, p) = jkj{p)fj{x, p), we have: 


Uj2{ipij{-,p)) ^ (Tji + crj2mj{X) 

Ujl{'lplj{-,P)) CFj 

where mj{X) = bj 2 {p)/bji{p) are the local Weyl functions on the rays TZj. Thus, (5.1) can be rewritten as 
follows: 

cTpi + ap2mp{X) gji + o-j2TOj(A) f/i2(^ir(-, p)) 

Now we note that all the Weyl functions j = 2,p — 1 are known since corresponding potentials qj were 
already recovered from the solution of problems IP{j),j = 2,p—l. Moreover, the Weyl-type solution 
il>ii{x,p) can also be found from the solution of the problem /P(l). Thus, all the terms in the right-hand 
side of (5.2) are known and we can use (5.2) to hnd the Weyl function nip that actually completes the 
solution of the problem IP{r). 

The following theorem summarizes our results. 

Theorem 5.1. Specification of the scattering data J determines uniquely the potential on T. The 
functions qk,k = l,p can be recovered by the following procedure. 

1. For k = l,p— 1 using the procedure 4.1 solve the problems IP{k) and recover qk- Calculate tpiiix^p) 
(while solving the problem IP{1)). 

2. Find nip from (5.2). 

3. Given rup recover qp{x),x > 0 by solving the (local) inverse spectral problem on the semi-axis cc > 0 by 
means of the procedure described in m- 
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